Shell-Filling Effect in the Entanglement Entropies of Spinful Fermions 
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We consider the von Neumann and Renyi entropies of the one dimensional quarter-filled Hubbard model. 
We observe that for periodic boundary conditions the entropies exhibit an unexpected dependence on system 
size: for L = 4 mod 8 the results are in agreement with expectations based on conformal field theory, while 
for L = mod 8 additional contributions arise. We explain this observation in terms of a shell-filling effect, 
and develop a conformal field theory approach to calculate the extra term in the entropies. Similar shell filling 
effects in entanglement entropies are expected to be present in higher dimensions and for other multicomponent 
systems. 
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Over the course of the last decade, entanglement mea- 
sures have developed into a powerful tool for analyzing many- 
particle quantum systems, in particular in relation to quantum 
criticality and topological order Qj. Within the realm of one 
dimensional (ID) systems, arguably the most important result 
concerns the universal behaviour in critical theories, which is 
characterized by the central charge of the underlying confor- 
mal field theory (CFT) J2J-3J. Let consider the ground state 
|GS) of a finite, periodic ID system of length L and parti- 
tion the latter into a finite block A of length I and its com- 
plement A. The density matrix of the entire system is then 
p = | GS) (GS|, and we will denote the reduced density matrix 
of block Aby pa- Widely used measures of entanglement are 
the Renyi entropies 
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They encode the full information on the spectrum of pa G), 
and in the limit n — > 1 reduce to the von Neumann entropy 
Si = — TvpA In pa- When the subsystem size I is large com- 
pared to the lattice spacing, 5„ are given by 
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where c is the central charge, c' n are non-universal additive 
constants, and o(l) denotes terms that vanish for £ — > oo. The 
result Q has been confirmed for many spin-chains and itiner- 
ant lattice models, see [ 1 1 for recent reviews. The knowledge 
of the entanglement entropies has led to a deeper understand- 
ing of numerical algorithms based on matrix product states [6 1 
and has aided the development of novel computational meth- 
ods 0. 

The Hubbard model is a central paradigm of strongly cor- 
related electron systems. Its ID version has attracted much at- 
tention for decades, because it is exactly solvable and exhibits 
a Mott metal to insulator transition |8 |. The Hamiltonian for 
periodic boundary conditions is 
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.Cj+i 1<7 +h.c. 
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where a are fermionic spin- 1 creation operators at site j 
with spin a =f, I, rij t(7 



Cj. a , and we will assume re- 



pulsive interactions U > 0. In the following we will for the 
sake of definiteness fix the band filling to be one electron per 
two sites, i.e. TY-j- = JVj. = ^, but we stress that our find- 
ings generalize to other fillings and, in fact, to other models. 
It is known from the exact solution that the ground state of 
Q below half filling (less than one fermion per site) is metal- 
lic and the low energy physics of the model is described by 
a spin and charge separated Luttinger liquid [ 8 1 equivalent to 
the semi-direct product of two c = 1 CFTs (9). 

Given this state of affairs, it is quite surprising that the en- 
tanglement entropies do not always follow This is shown 
in Fig. [T] which shows numerical results for 5i obtained by 
density matrix renormalization group (DMRG) for a quarter- 
filled Hubbard model at U = t for a number of different lat- 
tice lengths L. The data is seen to collapse on two curves, with 
only the lower one following the CFT result Interestingly, 
both the L = 4 mod 8 and the L = mod 8 data exhibit 
scaling collapse, but to different functions. We stress that this 
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FIG. 1 : DMRG data for Si - 2/3 In L as a function of x = £/L for 
U = t and L = 24, 28, 32, 36, 40, 44, 48, 52, 56, 60, 64. The lower 
and upper branches corresponds to lattice lengths L = 4 mod 8 and 
L = mod 8 respectively. 
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behaviour is very different from the lattice "parity effects" for 
Luttinger liquids ifTUl fTTl . which refer to o(l) corrections in 
S n > 2 only. 

The shell-filling effect. In order to understand the origin 
of the difference in entanglement entropies between L = 
4 mod 8 and L = mod 8, we consider the ground state 
in the limit U —> 0. Here we are dealing with with non- 
interacting, spinful fermions, for which the boundary condi- 
tions on a ring fix the momenta to be p m = 2-Kvn/L with inte- 
ger to € [— L/2, L/2). For a chain of length L = 8n+4, quar- 
ter filling corresponds to an odd number N a = L/A = 2n + 1 
of spin-cr fermion, and the unique ground state is the symmet- 
ric Fermi sea 



P = X)fLi kj> while the energy (in units of t) is given by 
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where c\{k) = L~ 1 / 2 Ylf=i e~ lkl Cj a are creation operators 
in momentum space and |0) is the fermionic vacuum state. 
On the other hand, when L = 8n, N a = L/4 = 2n is even 
and it is impossible for a given spin species to form a sym- 
metric Fermi sea. As a result the ground state is degener- 
ate. In particular, there are two degenerate ground states with 
N a = L/4 = 2n, that have zero momentum and are parity 
eigenstates (parity is a good quantum number) 



ct(fcp)C| (— kp) + ffcj (fcp)cl(— fcp) 
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Here kp = ir/4 is the Fermi momentum. As is shown below, 
the U — > limit of the Hubbard model ground state gives 
the state |+). The shell-filling effect is now clear: for L = 
4 mod 8 the ground state is a symmetrically filled Fermi sea, 
while for L = mod 8 it is given by the linear superposition 
of two asymmetric Fermi seas. In terms of spin symmetries 
this state corresponds to the S z = component of a S = 1 
multiplet. 

Be the Ansatz (BA) solution. We now turn to the case U > 0. 
Eigenstates of the Hubbard chain are parametrized in terms 
of the solutions { A Q , kj } of the following set of coupled BA 
equations ll8llT6l 
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Here u = U/{4t), 9{x) = 2arctan(x) and N = A t + N±. 
For real solutions of the BA equations (|6j, the "quantum num- 
bers" Ij (J a ) are integers if is even (if N+ is odd) and 
half-odd integers if N± is odd (if iV-f is even). The mo- 
mentum is expressed in terms of the parameters {A Q , kj} by 



E = uL 



N 
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where fi is the chemical potential. Following Ref. ifTTl . we 
define regular BA states as eigenstates of Eq. <|3j arising from 
solutions of <j6j with 2JVj, < N, where all kj and A a are finite. 
We denote these states by \{Ij}\ { J a })reg- As was shown in 
Ref. IfTTl . all regular BA states are lowest-weight states with 
respect to the SO(4) symmetry of the Hubbard model ifTHl . 
and a complete set of energy eigenstates is obtained by acting 
on them with the SO(4) raising operators. For L = 4 mod 8 it 
is known ll8l [T6l that the quarter filled ground state is a regular 
BA state characterized by the choice Ij — — 2n- 
1, . . . , An + 2 and J a — — n — 1 + a , a = 1, . 

For L = 8n (n a positive integer), we find that there are 
two degenerate lowest energy regular, real solutions of ([6]) 
with Nf = iVj, = 2n fermions. They are obtained by the 
two choices Ja' 2 ^ — — n — h + a , a — 1, . . . , 2n and 



t+3,3 
.,2n + l. 



-2n- 



(2) 



j + j ,3 = 1,.. .,4norPf = -2n-l+j , j - 

1, . . . , An. We stress that the distribution of the Ij is asym- 
metric around zero in both cases. Interestingly, these are not 
ground states. The regular solution with the lowest energy 
involves one pair of complex conjugate A Q 's known as a 2- 
string, but it is not the ground state either. 

Let us now consider regular BA states with total spin quan- 



tum number S z 



1, i.e. N- 



t 



2n + 1, N± = 2n - 1. 



These are by construction lowest weight states of the spin- 
SU(2) symmetry algebra. The lowest energy regular BA state 



in this sector corresponds to the (symmetric) choice P 
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1, . . . , 2n — 1. Crucially, the state 

S-|{/j 0) };{jW}) r e 
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is a (non-regular) eigenstate of the Hubbard Hamiltonian with 
Nf = N± = L/8 fermions. Here S~ = c ] i c j-t ^ s 

the spin lowering operator. As [S~ , H] = its energy is 
the same as that of the regular BA state { ji°'}) r cg- 

The energy difference between ([8]) and the regular solutions 
discussed above can be calculated for large L using standard 
methods [9| and is found to be negative. Considering other 
non-regular Bethe Ansatz states in an analogous way, we find 
that ([8]) is in fact the ground state. 

Bosonization. The low-energy physics of the Hubbard 
model is described by a spin-charge separated two-component 
Luttinger liquid Hamiltonian lfl2l 



(9) 



H = E y / dx i( d ^-) 2 + (^e a ) 2 ] , 



where v c , s are the velocities of the collective charge and spin 
degrees of freedom. For L = mod 8 the mode expansions 
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of the canonical Bose fields $ a = ip a 
fields 6 a = cp a — (p a follow from 



(p a , and their dual 



ip a (x,t) = P a + 



<p a (x,t) = P a + 
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where x± — x±vt and ao is the lattice spacing. The structure 
of the ground state for L = mod 8 is encoded in the zero 
modes, which have commutations relations [P a , Q a ] = — | = 
— [Pa., Qa]. The eigenvalues of Q a are 



9c = V8F E (i^c + l)m a + (1 - if c )r 
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where X c is the Luttinger parameter in the charge sector, mo- 
are half odd-integer numbers, and the eigenvalues of Qa are 
obtained by interchanging m a O m CT . The Hamiltonian then 
has the mode expansion 



H = Y ~T~ Ql + Ql + Y 2nn ( a hn a a,n+a{ n aa,n) 



(12) 
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There are two degenerate ground states 

l±> = ~[|1,0;0,1) ±|0,1;1,0)], (13) 

where we have introduced a notation \rrif, m^; m-j-, mj,) for 
states that are annihilated by all a a .„, a a , n and have eigen- 
values q a (m a ,rn a ) and qjjn a , fh a ) of the zero mode oper- 
ators Qa an d Qa. respectively. In the Hubbard model the de- 
generacy between |+) and |— ) is removed by the presence 
of a marginally irrelevant interaction of spin currents and the 
ground state in fact corresponds to |+). Carrying out the usual 
conformal map from the cylinder to the plane, we can express 
this state in radial quantization as lfT3l 



|+> oc lim cos(V27r$ a (z,£))|0) 



(14) 



where |0) is the vacuum state of the free boson theory ^ on 
the plane and z = exp {^-(vt — ix)), z = exp(-^-(vt + 
ix)). The key result of the above considerations is that in 
the Luttinger liquid approximation to the Hubbard model, the 
ground state for L = mod 8 is given by the excited state 



( 14 1 of the free boson theory on the plane. We note that the 



operator cos ^V27r$ s J is not a local operator on the cylinder. 

CFT approach to the Renyi entropies. The above consider- 
ations permit us to reduce the calculation of the entanglement 
entropies in the ground state of the Hubbard model to a CFT 
calculation of the entanglement in a low-lying excited state. A 
general approach to the latter problem has been recently de- 
veloped by Alcaraz et al. lfT4l[T5l and their main result can be 



summarized as follows. The n'th Renyi entropy for an excited 
state of the form O(0, 0) 1 0) is given by 
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where c' n is a O -independent constant, and the scaling func- 
tions F® (x) are given by 



F n {x) 



(m= O(^ + 2k))Ol(*-(-x + 2k))} 
n 2n ( h+h )(0(irx)Oi( - irx)) n 
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Here h and h are the conformal dimensions of the operator 
O. In our case 0(x) = 2 cos(^/2tt^ s (x)) and we need to 
evaluate (a = Ya=i a i) 



2n 



j — l CTi, ...,<7 n — ± l<j 
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where the x/s are given by ( 16 1. We find that F*(x) can be 



expressed as the square root of a determinant, which, surpris- 
ingly, is identical to Eq. (56) of Ref. [ 15 1. We have succeeded 
in expressing this determinant in a form amenable for analytic 
continuation in n 



[Fn(x)} 



71 



(n-2p+l) 2 



sin 2 (7rx) 



P =i 



2 sin(7ra;) 



l+n+n csc(7ra:) 
2 



I l — n+n csc(ttx)\ 

V 2 J 



■ (18) 



Using that for the Hubbard model c — 2we then obtain a CFT 
prediction for the shell filling effect by combining equations 



( 18 1 and (15 1. In order to obtain an expression for the von 
Neumann entropy we need to take the limit n — > 1, which 

gives 



c (L irl 
bx = - In — sin 

3 V 7T 
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+ c[ - F{(x) + o(L) 
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where 

F[{x) = ln|2sin(7ra;) 



1 

2 sin(7ra;) 



+ sin(7T2;). (20) 



Here ip(x) is the digamma function. We note that both ( 18 1 



and ( 19 1 apply also to certain excited states in spin chains |[T5l . 
For small x, we have (Fi(x))' — ir 2 x 2 /3 + 0(x 4 ) in agree- 
ment with the general result in 1 14 1. 

Comparison with numerical results. We performed exten- 
sive DMRG lfl9l computations of the periodic quarter-filled 
Hubbard model by keeping M — 3000 states in order to 
achieve satisfactory convergence for periodic systems up to 
length L = 64. For small values of U < i we find good 



agreement for both Si and 52 with the predictions (19i and 
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FIG. 2: SSi = Sx - | In [£ sin - cj as a function of x = 

^/i for (7 = 0.3t and L = 24,32,40,48,56,64. The constant 
c[ — 1.205 has been adjusted by hand. The solid curve is — F{(x), 

(jT3J. A representative example is shown in Fig. [2] As ex- 
pected the agreement with the CFT prediction is best for large 
block lengths £ ~ L/2 and becomes poor for small t, when 
lattice effects become important. In this region Si further- 
more exhibits strong oscillatory behaviour as expected iflOl . 
For larger values of U > t the agreement with the CFT pre- 
dictions for both Si and S% becomes increasingly poor. We 
now turn to the origin of these discrepancies. 

Effects of the marginal perturbation. It is well know that in 
the Hubbard model the low-energy Luttinger liquid Hamilto- 
nian Q is perturbed by a marginally irrelevant operator in the 
spin sector 11121 . This leads to logarithmic corrections ||20l , 
which can be quite important for small system sizes. The ef- 
fects of a marginal perturbation on the ground state entangle- 
ment in CFTs was studied in lETl . These corrections are small 
for the isotropic Heisenberg chain 11221 as well as the Hubbard 
model for L — 4 mod 8. However, the effects of the marginal 
perturbation on the shell-filling effect are quite large already 
for moderate values of U > 2t. In order to quantify them, we 
have considered an extended Hubbard model 

H ext = Hn uhh + V 2 ^ n i,o n j+2,o'- ( 21 ) 

At weak coupling the main effect of V2 is to reduce the bare 
coupling constant of the marginal perturbation. We note that a 
nearest-neighbour density-density interaction would be inef- 
fective at quarter filling and weak coupling 11231 . We find that 
increasing V2 from zero leads to a significant improvement in 
the agreement between the CFT prediction (TT9J, (JT3J for the 
shell-filling effect for the available system sizes L < 64. 

Conclusions. We have described a novel shell-filling effect 
in entanglement entropies of the ID quarter-filled Hubbard 
model with periodic boundary conditions. We have developed 
a CFT approach to calculate the additional contribution to the 
Renyi entropies, and found good agreement with numerical 
computations. The effect, while somewhat unexpected, has 
a simple origin: for certain ratios of lattice lengths to par- 
ticle numbers in multi-component systems, the ground state 



cannot be thought of in terms of a product of Fermi seas (in 
general these will consist of appropriate elementary excita- 
tions), but is in fact a linear combination of different such seas. 
This suggests similarities with the results obtained [24 25 1 for 
the entanglement of linear combinations of degenerate ground 
states. However, in our case the ground state is unique for 
U > (and fixed S z = 0) and is thus not based on a degener- 
acy. Finally, we expect shell-filling effects to exist for interact- 
ing bosons as well as for fermions in one dimension, as well as 
in higher dimensional critical systems. They can for example 
play a role in numerical studies of two-dimensional gapless 
spin liquids, which display a spinon Fermi surface 126H281 . 
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